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Abstract — A comprehensive theory for robust PID control in 
continuous-time and discrete-time domain is reviewed in this 
paper. For a given finite set of linear time invariant plants, 
algorithms for fast computation of robustly stabilizing regions 
in the (kp , ki, fco)-parameter space are introduced. The main 
impetus is given by the fact that non-convex stable regions in 
the PID parameter space can be built up by convex polygonal 
slices. A simple and an elegant theory evolved in the last few 
years up to a quite mature level. 

I. INTRODUCTION 

It is a well-known fact that by far the most applied 
control law for SISO systems in nearly all industries (process 
control, motion control, aerospace etc) is the PID control. 
No other controller can compete to PID when it comes 
to performance per simplicity of the control structure, this 
being the reason for its absolute dominance in the practice. 
Traditionally, tuning of a PID controller has been the over- 
whelming usability approach in research and applications. 
The design technique presented here is in some sense an 
opposing one. We want to compute the total set of PID- 
stabilizers. While it turns out to be interesting in theoretical 
terms, its impact on practical applications is difficult to 
predict. Advanced software tools based on this technology 
(e.g. ROBSIN, [Bajcinca and Hulin (2004)]) suggest to the 
user a 3D-region in (kp, k], &£>)-parameter space, where 
he can select a controller from. By doing so, he would 
additionally have an idea how robust (i.e. how far from 
instability) his design is. A further good news is that the same 
technique applies for time-delay systems and in discrete-time 
domain. 

This paper covers theoretical design issues only. It has 
been shown that the stabilizing region for continuous- 
time PID (= kp + k] /s + kp>s) controllers is defined by 
a set of convex polygonal slices normal to kp axis in 
the (kp, kj, kp,) parameter space, see [Ho et.al. (2000)], 
[Ackermann and Kaesbauer (2001)], [Bajcinca (2006)]. The 
method followed by the author uses the decoupling ef- 
fect of PID parameter space at singular frequencies, 
[Bajcinca (2006)]. Thereby the characteristic polynomial de- 
couples into two frequency parameterized equations, one 
involving ki and kp,, and another one with kp only. As a 
consequence, the computation of all stable PID controllers 
may be divided into two subproblems: assertion of stable in- 
tervals of parameter kp (so-called £p-problem), and detection 
of stable polygons on the plane (kj, kp,) for a given kp. In 
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[Ackermann et.al. (2002)] and [Bajcinca (2007)] the design 
approach was transferred to discrete-time systems, and in 
[Bajcinca (2004)] to time-delay systems. 

A general algorithm for automatic detection of stable poly- 
gons is based on the analysis of the motion of eigenvalues 
in the vicinity of the singular frequencies. This algorithm 
was first presented in [Bajcinca (2001)]. It detects the inner 
polygons and selects the one with maximal stable eigenval- 
ues, which is finally checked for stability. It is important to 
emphasize, that this algorithm can be applied equally well 
for PID quasi-polynomials describing the PID control loops 
for time-delay systems. 

Preserving simplicity has been a basic motivation in 
searching for a solution to the £/>-problem. A simple neces- 
sary condition was firstly understood in [Bajcinca (2005)]. It 
turns out that for a given plant a sufficient number of singular 
frequencies must be available for its stabilization. Since kp 
uniquely determines the number of singular frequencies, one 
can directly read ^-intervals from the kp-plot, which may 
host stable polygons. This simple criterion can be extended 
to discrete-time domain and time-delay quasi-polynomials, 
for instance, see [Bajcinca (2007)] and [Bajcinca (2004)]. 

The presented methods may be directly applied for the 
computation of the total robust set of PID controllers which 
stabilize a finite set of plants (multi-model). This is exactly 
what is meant by robust design of PID controllers. 

In this paper methods for PID control in continuous- and 
discrete-time domain, as well as for time-delay systems, 
are reviewed. Therefore, we will have to switch between 
the s domain (continuous-time) and z-domain (discrete-time) 
while presenting results. The paper is organized as fol- 
lows. In Section we postulate the stabilizability problem. 
Section 



III 



presents all steps of the design process, and 
simultaneously introduces the formulates technical problems, 
in a simple-case study. Theoretical fundamentals of the 
methods are introduced in Section [IV] Automatic detection 
of convex stable polygons and k — P -problems are discussed 
in Section [V] and VI respectively. Finally, in Section VII the 
ideas are extended for time-delay systems. For illustration 
purposes, throughout the paper examples are provided. The 
reader may follow them by using the toolbox ROBSIN, which 
can be downloaded from http://www.robotic.dlr.de/robsin. 

This article primarily reviews the original work of the 
author. It has not been an intention to refer to the all work in 
the area. Still, the key contributions (and contributors!) to the 
theory are referred to and they read as follows. The role of 
convex polygons was firstly understood by Battacharyya and 
his co-workers, see the monograph [Ho et.al. (2000)]. Their 
derivation bases on Hermite-Biehler theorem. For computa- 



tion of stable polygons they proposed linear programming 
techniques, but they did not really address the ^-problem. 
Munro's computation is based on the real-axis intersections 
of the Nyquist plot, see, [Munro and Soylemez (2000)], 
[Soylemez, Munro and Baki (2003)]. The relationship to 
singular frequencies was noticed firstly by Ackermann and 
Kaesbauer, see [Ackermann and Kaesbauer (2001)]. Soyle- 
mez proposed a solution to the £p-problem, however in 
author's opinion his criterion is not as simple and usable 
as the one presented in this article. The remainder theory 
reviewed here is fairly developed by the author, including 
the algorithms for automatic detection of stable polygons 
and solutions to the £p-problem, as well as generalizations 
to discrete-time domain and time-delay systems. 

II. Problem definition 

A. Discrete-time domain 

Consider a simple closed curve dF = {z \ z = t(oc) + 
j T}(a), a is a parameter}, in z— plane, which is symmetric 
to the real axis T and can be expressed in the form 



Let 



ST: F(t,tj)=0. 



p=A( z )Q(z,n,r 2 ,r 3 )+B(z)=0 



(1) 



(2) 



be a three-term algebraic equation in (r\, r2, r 3 ), with given 
polynomials A and B 



A(z) = ao+flizH Ya m i n 

B{z) = b +b lZ +---+b n z n , 

and the second-order polynomial Q in the form 

Q = Si(z)n + 8 2 {z)r 2 + 8 3 {z)r 3 . 



(3) 
(4) 



(5) 



In this article we want to compute the set of all parameters 
r\, r 2 , J"3, s.t. the polynomial Q is T— stable, that is, all its 
eigenvalues must lie within the T-region (enclosed by dF). 
Of main interest are circles with center on the real axis T and 
an arbitrary radius, which will be referred to as T— circles. 
For discrete-time systems particularly important is the unity 
Schur-circle. 

It may be easily shown that |2|i describes the characteristic 
equation of a feedback loop with a PID or a three-term 
controller. Indeed, a discrete-time equivalent of the PID 
controller has the transfer function 



C\ +C2Z + C3Z 

(z + zi)(z-l) 



(6) 



Its structure follows in the quasi-continuous consideration by 
applying the rectangular integration rule (s— > (z— \ )/Tz) to 
the ideal PID controller kj/s + kp + kps, resulting in z\ =0, 
or by the trapezoidal integration rule (s —> 2 (z — 1 ) /T (z + 1 ) ) , 
resulting in zi = 1. Also the realizable PID controller ki/s + 
kp + kos j (1 + T[s) converts by the trapezoidal integration 
rule to the controller structure (|6]i with a pole at z\ — 



— (2T\—T)/{2T\+T). Equation (|2j covers also a three- 
term controller with an arbitrary second order denominator 
polynomial 

*(z) = " (z)(Cl + C f W) . (7) 
d(z) 

For both control structures (|6]l and (|7]), the polynomial Q is 
of the form 

Q = c l +c 2 z + c 3 z 2 . (8) 

It is clear that |5]) and (|8]l are connected by a linear trans- 
formation 

c = T r, det T ^ 0, (9) 

with c = [ci,C2,c 3 ] T , r = [ri,r2,r 3 ] T . The matrix T is 
determined by the polynomials d\(z) 1 o\(z), S 3 (z). 

B. Continuous-time domain 

The pendant equation |2]) for a feedback-loop with a PID 
controller in continuous-time domain is 



p =A(s)(k I + k P s + k D s 2 )+B(s) 



(10) 



where polynomials A(s) and B(s) are as in <j3j and Q. As 
in the discrete-time case, we want to compute the set of 
all parameters kp, kj and kp> for which the polynomial ( |T0] > 
is Hurwitz-stable, i.e. all of its eigenvalues lie on the left- 
hand side of the imaginary axis s = j(0. In other words, 
here <9r = {jco : Vo) e R}. In view of the definition the 
polynomial Q is also of the form ([8} 



Q = kj + kps + kp>s 



(11) 



Obviously discrete-time consideration is more general, with 
z,(X,ri,r2,r 3 corresponding to s,co,ki,ko, kp, respectively. 

III. The very basic idea 
Consider the special case with A(s) = 1 in ( 10 1. Substitu- 



tion s = j(0, decouples the latter into two equations 

ki-0) 2 k D = -R B , k P = -I B /(0 



(12) 



where R and / stand for the real and imaginary part of 
polynomial B(s) at s = jco. Notice that PID parameters 
appear decoupled in tow equations. Computation of the stable 
PID region should be quite obvious for this case study. (1): 
First, for a fixed parameter kp = k' P solve for the frequencies 
CO 1 from the second equation in (12i, representing the kp- 



plot and shown in Fig. 1(a) Such frequencies are known as 
singular frequencies. (2): Map all singular frequencies using 



the first equation in (12i into the kp,)— plane as shown 



in Fig. |l(b)| and compute the stable polygons (gray area in 
Fig. [I}. Hereby, each pair of singular frequencies ± jco' maps 
to a straight boundary line. Note that, as parameters kj and 
kp> are varied, with kp kept fixed, the eigenvalues of Q can 
cross over imaginary axis jco at singular frequencies only. 
This procedure is repeated for other kp parameters, which 
yield other stable polygons. Thus a tomographic 3-D picture 
results, as that in Fig. [7] 

In the sequel, we want to generalize the above algorithm 
and provide solutions to the problems: [PI] For what values 
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(b) Mapping of singular frequencies into the (k[ , ^o)-plane 



Fig. 1. Illustrating the basic algorithm for computation of stable PID 
controllers in continuous-time domain 



of the parameter kp stable polygons should be searched for - 
indeed one would like to exclude a priori ^-intervals with 
no stable polygons. This problem is referred to as the kp- 
problem. [P2] For a given set of boundary lines, how to 
automate the computation of the stable polygons. 

IV. Basic definitions and theorems 

Let H and G represent the real and imaginary part of the 
characteristic polynomial p(z,ri,r 2 ,r 3 ) in Q on <9r. 

Definition 1: Y is said to be singular with respect to 
parameters r\ and r2 in (j2| ;/ 

d(H,G) 



Rank- 



1 for all z € <3r. 



(13) 



d{r\,r 2 ) 

The latter equation is referred to as the rank-condition. 
From now on, we only consider Ts which satisfy (13i. A 



zero of the polynomial Q that fulfills ( 13 1 is referred to as 
singular (or critical) frequency. 

Definition 2: A function E^(z) defined as 



with 



Q(z,n,r 2 ,r 3 ) = E r (z) q(z,n,r 2 ,r 3 ) (14) 
I q = r 3 g 3 (a)+go(a), ae[a,b] (15) 



where I q stands for the imaginary part of q, will be referred 
to as decoupling function of Q on dT. 

In other words, by introducing Ej-, a function q is extracted 
from Q, whose imaginary part depends on one parameter r 3 



only. It can be easily checked that if F satisfies ( 13 1, then 

dl a 



n dI 4 

dr\ dr 2 



o, Vz g <?r. 



(16) 



Furthermore, two trivial decoupling functions on <9r of Q in 
Q are 

E T (z) = 8i(z), E T {z) = &i{z). (17) 



Using these two facts, the next statement follows directly. 
Theorem 1: Consider the function 

P(z) 



F(z) 



(18) 



A(z)E T (zY 

The equation F(z) = for z G <3r decouples the parameters 
r\ , r 2 and r 3 into two equations 

rihi(a)+r 2 h 2 {a)+h Q {a) = 0, (19) 

r 3 g 3 (a)+go{a) = 0. (20) 
Note that the latter equations reveal decoupling of the 
parameter space (ri,r 2 ,r 3 ). They represent generalizations 



of the simple equations ( 12 1 



Finally, let {p v } represent a finite set of polynomials of 
the form |2]), i.e. 



Pv =A v (z)Q{z,ri,r 2 ,r 3 )+B v (z). 



(21) 



For instance, this may be a multi-model of a continuum of 
plants or Kharitonov polynomials of an interval uncertainty. 
It can be simply proven that the rank-condition ( fT3"] l does 
not depend on the polynomials A v (z) and B v (z). Hence, a 
singular T is completely determined by the polynomial Q in 
(15). The polynomials A v (z) and B v (z) rather determine the 
singular frequencies on <9r. 

A. Hurwitz-stability 

Consider the Hurwitz-region dT 
condition (fT3]l 

d(H,G) 



Rank-. 



--{jo: co G R}. Then 
1 (22) 



d{k h k D ) 

is satisfied everywhere on the imaginary axis, since 



H = R A ki-(O z R A k D -(Ql A kp+R B (23) 
G = I A ki-co 2 I A k D + coR A kp+I B . (24) 

Referring to equations ( [T7j ) and (111, a simple choice for the 
decoupling function is Er(z) = 1. Indeed the function ( 18 1 

A(s) 

on the imaginary axis s — jco decouples into two equations 



F(s) — ki + kps + kos 1 



(25) 



of the form ( 19 1 and ( 20 1, similar (but not equal) to equations 
( [12"} , Note that the rank-condition p2| ) applies also on any 
line parallel a — Oq to s — j(0. Thus all derivations hold also 
for <?r = {fj + j(0 : co G K}, see [Bajcinca (2006)]. 

B. Schur-stability 

Consider the Schur-circle <9ri = {e' a : a G [-7C, n]}. It 



can be easily checked that the rank-condition ( 13 1 is satisfied 
on dF\ for 

Q = (l+z?)n+zr 2 + r 3 
and that the matrix T, as defined in |9|, reads 

1 1 

1 

1 



(26) 



(27) 




Fig. 2. Schur- and T— circles 



Following ( 17 i, the trivial decoupling functions of (26i on 
the Schur-circle F\ are 

E r {z) = l+z 2 or E r (z)=z. (28) 

For Et(z) = z, the imaginary part of the function on Tj 

F(z) = — n+r 2 + ^r 3 (29) 
z zA(z) 

is of the form ( jljj l, since imaginary part of is null on 
C. Y-stability 

Consider a T— circle with center on real axis dT = {m + 
pe ja , a G [-71, %]}, Fig. [2] Now 



Q= (p 2 - m 2 + z 2 )ri + (z- m)r 2 + r 3 . 



(30) 



For T— circles with center at T = ra and radius p a transfor- 
mation matrix from c— to r— parameter space is 



p--m z 

1 



— m 1 
1 




Furhtermore 



Er(z) =z — m, 



and 



p 2 -m 2 +z 2 

F(zJ = n + r 2 - 

z — m 



B(z) 



(z — m)A(z 



(31) 



(32) 



(33) 



V. Stable convex polygons 

This section recalls briefly a solution to problem [P2] as 
formulated in Section [ill] For details the reader is referred 
to [Bajcinca (2001)] and [Bajcinca (2006)]. The algorithm 
is motivated by the concept of inner polygons, which claim 
a necessary condition for stability: A polygon 11 is said 
to be an inner polygon if entering the polygon in (r\,r 2 )- 
i.e. (fc/,fco)-plane, causes an eigenvalue-pair to enter the 
T— region at the corresponding singular frequency. 



A. Discrete-time domain 

Let {z 1 } be the set of singular frequencies on dF deter- 
mined by the equation ( 20 1 for a fixed r$, and let { X — X (z!) } 



be the corresponding straight lines determined by ( [T9| , see 
Fig. [3] In order to automate the detection of an inner polygon 
each boundary line A, will be assigned a "transition" e: it 
is negative if the transition [8r\,8r 2 ] (see Fig. [3]) over the 
singular line causes an eigenvalue to become stable, i.e. enter 
the T-region, otherwise it is positive. Let e\ correspond to 
8ri > 0, 8r 2 = and e 2 to Sr 2 > 0, 8r\ = 0. 



r 2 





Fig. 3. Definition of e\ and e^. the motion of eigenvalues in the vicinity 
of a singular frequency z'. The shaded part refers to the stable side of Y 
and the normal N points outside the T-region. 



To compute the functions e\ and e 2 , define the normal 
vector N on dT at the singular frequency z' = t(oc') ± jr\ (a') 
by its complex associate as 



N={dF/dx + jdF/drj), 



(34) 



For tracking the motion of an eigenvalue due to small 
variations in r\ and r 2 introduce the vectors 



Ml = 



dz 
dri 



dT 



■ dr l A 

■j— and ju 2 

dri 



d^ 
dr 2 



Assuming dp/dz ^ it can be easily shown that 
Mi 

Now transitions functions can be computed by 



dp ,dp 
dr\ dz 



h/2 = Re (MiVv) ■ 



(35) 



(36) 



(37) 



Using this information, an algorithm for the detection of 
the inner polygons (those with maximal number of T- 
stable eigenvalues) is developed in [Bajcinca (2001)]. Such 
polygons are the only stability candidates, that can be proved 
by checking any point therein. 

Example 1: Consider the discrete-time model of the plant 



10" 



5 4. 1 65z 3 + 45 ,77z 2 + 45 ,77z + 4. 165 
z 4 - 3.985z 3 + 5.97z 2 - 3.985z+ 1 

and a three-term stabilizer 



C(z) = 10' 



(z 2 - 1. 54k + 0.5992) (q +c 2 z + c 3 z 2 ) 



z(z + 0.4047) (z + 0.2162) (z - 0.4934) 
'For situations with dp/dz = 0, refer to [Bajcinca (2006)]. 



(38) 



(39) 



whose parameters ci,C2,C3 are to be synthesized. The synthesis is 
done in (r\ , ri, ^-parameter space. Therefore the transformation 
(27} can be used. Then 



0.59 



A(z) = z 5 +9.44z 4 -5.34z 3 -9.34z 2 + 5.04z-^ 
B(z) = 0.19z 8 -0.73z 7 +z 6 -0.45z 5 -0.12z 4 
0.14z 3 -0.009z 2 -0.008z. 



For n = —0.26118, the singular frequencies lying on the Schur- 
circle are computed to be z'j = 1,4 = 0.9172 ± j0.3983,z' 3 = 
0.5628 ± jO. 8266, z' 4 = — 1. The resulting stable polygon is shown 
in Fig. [5] It is enclosed by the straight lines frequencies X\ , A2 and 
A3, corresponding to z'pZ^ and z' 3 - 



''2 





stable polygon 















0.6 0.8 1 1.2 1.4 1.6 1.8 2 ; 

Fig. 4. The stable polygon for r$ = —0.26118, Example 1 



B. Hurwitz conditions 

Let ei correspond to 8ki > 0, 8kp> = and ep> to 8ko > 
0, 8ki = 0. It can be shown that (137} yields 



ei/D 



d(H,G) 



d(co,k I/D ) 



(40) 



Expressions for ej/ D do not depend on where a boundary 
line is crossed at. Indeed, check that (ftp} is equivalent to 



ei/D : 



dk P 
dco 



d(H,G) 



d(k P ,ki/ D ) 



(41) 



where kp = kp(co) stands for the fep-plot. Since co' is a 
singular frequency, the determinant in the latter equation is 
shown to be independent on kj and kpj- Furthermore, the 
following holds 



sign ei 



-sign 



dk P 
dco 



sign e D = sign 



dk P 
dco 



(42) 



These expressions prove again that the transitions e/ and ep> 
are independent on parameters kj and kp>, and of opposite 
sign. Their sign is determined by the slope of the fcp-plot at 



the corresponding singular frequency, see Fig. 1(a) 
Example 2: Consider the polynomial \\()\ with 

A(s) = -0.5.? 4 -7.? 3 -2i'+l 

B(s) = 5 7 + ll5 6 +465 5 +95/ + 109.s 3 +74i' 2 + 24.s. 
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Fig. 5. Stable polygons for kp = —2, Example 2 



The singular frequencies for kp = —2 are computed from its kp- 
plot (see Fig. [6j: = 0, s[ = ± ;0.3530, 4 = ± j0.6638, s' 3 = 
±y'0.7742, = ±y'3.3473. Fig.|5]depicts the corresponding straight 
lines and the stable polygons. Note that the stable polygons need 
not to be connected. 

VI. ^/.-Problem 

A. Hurwitz-stability 

This section focuses on the problem [PI], as defined in 



Section III a rule is sought to discriminate £p— intervals 
with stable PID controllers. Intuitively, it must be tightly 
related to the £/>-plot. Indeed, it is clear that at maxima 
and minima of the kp— plot - compare Fig. 1(a) and 1(b) 
- convex polygons close at an edge due to the overlapping 
of two straight boundary lines. As ^-intervals with maximal 
number of singular frequencies are most likely to host stable 
polygons. The following result renders this idea precisely. 

Theorem 2: Consider the polynomial \10\ . Assume that 
polynomial A{s) has no zeros o n th e imaginary axis and let 

N: order of the polynomial ( 10 \ 

M: order of the polynomial A(s) 

P: number of RHP zeros of A(s) 

Z: number of singular frequencies in the interval 

< co < +°°. 
A necessary condition for stability of ( |70| > is 

z ^E(N-M + 2P-l)_ (43) 

Here, the function £:Nh^ selects the nearest smaller 
even number. The proof of the theorem can be found in 
[Bajcinca (2006)]. It is important to observe, that we exclude 
the zero singular frequency from Z. 

Using this criterion one can directly read from the kp-plot 



(Fig. 1(a) i the kp— interval(s) with (potentially) host stable 
polygons. However, in some cases a polygon may close when 
three boundary lines intersect at a single point in (kp,k],kpj)- 
parameter space. This situation destroys the sufficiency of the 
condition (see Lemma [T} is not sensed by the above criterion 
and will be discussed in Section IVI-CI 

The following is an extension of Theorem [2] to the cases 
when A(s) possesses zeros on the imaginary axis. 



10 



CO, 



-25 



Fig. 6. The i/>-plot, Example 2 



Theorem 3: Suppose A(s) has J zeros on the imaginary 



axis. Then, for stability of the polynomial (701 Z, singular 
frequencies are required within the interval CO € (0, +°°), 
where 

(a) if s — is not a zero of A{s) 



Z> 



E{N-M + 2P-J- 1) 



(b) if s = is a zero of order Jq of A(s) 

z> E(N-M + 2P-J-1)-E{J ) 



(44) 



(45) 



stable polygonal slice can be computed. The result is shown in 

Example 3: Separated stable kp— intervals. Consider the poly- 
nomial < | 1 0| > with 



A(s) 
B(s) 



s 3 + 3s 2 - 
5 5 + 2, 4 - 



-3s J +7s z + 14s. 



It can be directly read that N = 5, M = 3 and P = 2. According to 
Theorem [2] for stability at least 2 positive singular frequencies are 
required. Now consider Fig. [8] where the kp-p\ot for CO > is de- 
picted. Two fcp-intervals of interest are directly recognized, namely 
— 1.8708 < kp < —1.5556 with 3 positive singular frequencies, and 
0.3157 <kp < 0.5333 with 4 ones. For other kps no stable polygons 
exist. 




Fig. 8. The kp-p\ot, Example 3 




15 



Fig. 7. The region of PID stabilizers, Example 2 

Example 2: (cont). The kp-p\ot is depicted in Fig. [6] Notice that 
it is very easy to read the number of available singular frequencies 
from Fig. [6] for a given kp. For the polynomial A(s) we have 
N = 7, M = 4, P = 1, J = 0, J = 0. According to Theorem [2] 
for stability at least Z > E(N - M + 2P - J - 1) /2 = 2 singular 
frequencies are required for ft) > 0. By checking the Fig. [6] it is 
obviously that this condition is fulfilled for —24 < kp < 6.1565. 
More precisely, for -24 <k P < -2.7614 and 3.7664 <k P < 6.1565, 
two non-zero singular frequencies exist, and for —2.7614 < kp < 
3.7664 four ones. Finally, by gridding kp within these intervals, 



Example 4: Missing stability. Let 
A(s) = 1 

B(s) = s 5 +s 4 -3s 3 -s 2 + 2s. 

Theorem [2] requires at least 2 singular frequencies in CO > 0, how- 
ever for —2 > kp, 1 singular frequency exists, otherwise none. Thus, 
polynomial \10\ is always unstable, no matter what kp, kj, kr>. 

B. Schur-stability 

Without loss of generality, we consider just the Schur- 
circle. The generalizations for other T— circles are straight- 
forward. 

Theorem 4: Consider the characteristic polynomial ([2]) 
and the Schur-circle T\. Let 
N: order of the polynomial (|2j) 
R: number of zeros of A(z)Er(z) lying inside dT\ 
J: number of zeros ^ ±1 of A(z)Ey{z) lying on dY\ 
J + : order of the zero +1 of A (z)Er(z) 
J : order of the zero —1 of A (z)Er(z) 
Z: number of singular frequencies in the interval 

< a < +7t. 
A necessary condition for stability of (J3J is 

Z>N-R- J + E{J+)+ 2 E{J - )+1 . (46) 

Example 1: (cont). It can be checked that A(z) has three zeros 
inside the Schur-circle, one zero at z = — 1 and one zero outside 
the Schur-circle. Thus if the decoupling function Er(z) = z is used, 



it follows that N = 8, R = 3 + 1 = 4, 7=1, /+ = 0, and /_ = 1. 
Hence, for stability, Z > 3 singular frequencies are required in the 
interval < a < +71. In order to discriminate stable intervals 
check the fcp-plot in Fig. [5] The stable interval is indicated by 
the grayed strip in Fig. [9] -0.52236 < r 3 < 0.00290. Notice that 
the zoomed plot in Fig. 15] reveals that for < rj < 0.00290 four 
additional singular frequencies appear. 

On the other side if the decoupling function Er(z) = 1 + z 2 is 
used, then N = 8, R = 3, / = 3, J+ = 0, and /_ = 1 i.e. again 
for stability Z > 3 singular frequencies are required in the interval 
< a < +7i. 




1 DC 15 



Fig. 9. The kp-plot. Example 1 



Example 5: PID control. Consider PID control of the plant {38} 
now using the control law j6|. It can be shown that in this case 



A(z) 
B(z) 



z 3 + 10.98 z 2 - 



10.98z+l 



O.^-O^+z 4 



-z 3 +0.5z 2 - 



O.lz. 



(47) 
(48) 



By using £r(z) = z, it is easily checked that N = 6, R = l + l = 
2, J = 0, J + = and /_ = 1. Hence, Z > 3 singular frequencies 
within < a < +7T are required. However, the maximal number of 
singular frequencies within < a < +71 is 2, so no PID controller 
can stabilize the plant l |38[ l. 

C. Stability peaks 

Example 6: Stability peaks. Consider polynomials 



A(s) 
B(s) 



1890/ + 658i+215 



1032 



s' + 



25 

92785263 4 



6175327 fi 98620159 , 

•s + r + 



r + 



10000 
97588159 3 



r + 



25000 
5413746 2 



10000 25000 625 

According to definitions in Theorem [2)^ = 8, M — 2, P = 0. 
Therefore for stabilizability Z>E(N — M + 2P— 1 )/2 = 3 singular 
frequencies are required in the interval < ft) < +°°. It can be 
shown that for both, hp = —9 and kp = — 10, Z = 3 positive singular 



frequencies are available, which satisfy the stabilizability condition. 
But, for kp — — 9 a stable polygon exists, and for kp = — 10 it does 
not. In other words, for some kp in between, the stable polygons 
must close in a vertex (and not an edge). To handle this situation 
one has to detect stability peaks within the intervals provided by 
Theorem [2] 

Note that at a peak, the three-term polynomial (JTOj must 
have at least three different eigenvalues on the imaginary 
axis, that is 



A 0) (ki + k P s + k D s 2 ) + B (s) 



■ R(s) J] (* 2 

(=1,2,3 



col 2 ), (49) 



where R(s) is a remainder polynomial, which has to be 
stable, otherwise the peak is irrelevant. This is a nonlinear 
set of ./V equations with N unknowns. The left-hand side of 
the equation provides the three unknowns kj,kp,ko and the 
right-hand side the rest N — 3 ones, including three singular 
frequencies 03- , i = 1,2,3 and N — 6 coefficients of the 
polynomial R(s). Hence by elimination of the latter N — 3 
variables, a system of three nonlinear equations with the three 
kj,kp,ko variables results. Its solution provides the required 
peaks. In general, finitely many stability peaks exist. 

Example 6: (cont). It can be shown that for the three-term 
polynomial defined by \A9\ a stability peak appears at kp ~ 
—9.0023. The three straight lines corresponding to <a[ ~ 0.2581, 
a' 2 w 0.44261 and co^ w 9.7621, intersect at k D w 21.4958, h w 
3.0195. A stability peak appears also in Example 7, see Fig. |10| 

Lemma 1: The condition in Theorem [2] is also sufficient 
for N < 6. Theorem [2] provides necessary and sufficient 
conditions for any PID feedback loop with plants of 1st, 
2nd and 3rd order. 

VII. PID FOR TIME-DELAY SYSTEMS 

In this section we extend the PID control theory to systems 
with time-delay. Such feedback loops always result in quasi- 
polynomials of the form 

2\ I D/^N^Xs 



p = A(s)(ki + k P s + k D s ) +B(s)e 



(50) 



where L > is the time-delay. Fundamental stability 
conditions regarding quasi-polynomials are provided in 
[Pontryagin (1955)]. E.g. a simple necessary condition is 
the existence of principal term, i.e. e in (50 1 must be 
multiplied by the highest power in s. Thus, in this section we 
assume n > m + 2, i.e. the quasi-polynomials of the retarded 
(n > m + 2) and neutral type (n = m + 2) are considered only. 
It is easy to check that decoupling conditions hold for the 



quasi-polynomial (50 1, too. With definitions 



a (fl)) = 



IRb 2 +Ib 2 
Ra 2 +Ia 2 " 



tan0(ft>) = 



Ra h — I a Rb 
Ra Rb + IaIb 



equations (19i and (20 1 take the form 

ki-co 2 k D = a(co)cos(coL + (j)(co)), 
cokp = a((D)sin((oL + 0(a))). 



(51) 



(52) 
(53) 



Obviously kp-plot is now a sinusoidal function and we have 
to deal with infinitely many singular frequencies. 



A. High-frequency behavior 



Consider equation ( |5T[ ). Given that n > m + 2 for high 

(54) 



singular frequencies, i.e. as ft) 



alco) 



co 



and 



tan ^ (ft)) — » 



±c 





0(ffl) 



or 71. 



For a fixed £p-grid, equation ( |53| ) implies 

1 



sin(ft)L + 0(fi))) 



-ft) 



l-n-l 



0, 



ft) 



that is, all higher singular frequencies tend to 
kn/L, 

{k+\)n/L, fceN,&>l, 
and 

cos(ft)L + </>(&))) -» ±1. 



(55) 



(56) 



(57) 



(58) 



To investigate the behavior of boundary lines for high 
frequencies, applying ft) — » °° to ( 52 1 reads 



&d = ±a(ft))/ft) 2 



(59) 



For quasi-polynomials of neutral type, n =m+2 in ( |54| >, and 
boundary lines converge to so-called infinity root boundaries 



k D = ±b n /a„ 



(60) 



For quasi-polynomials of retarded type, straight lines diverge. 
Note that infinity root boundaries describe the situation with 
infinitely many eigenvalues arbitrarily close to the imaginary 
axis. 

B. Relevant frequency range 



According to (42 1 the sign of the transition function e^ D 
at a singular frequency ft)' is determined by the slope of 
the function kp = kp(co) at ft)'. Hence, at successive singular 
frequencies corresponding to a fixed kp, e// D takes opposite 
signs. This motivates the definition of the set of odd £l and 
even Q. e singular frequencies 

• }, Q e = {<< 04, ■••} (6i) 



Q. = {CO[, CO3, CO5, 
with 



= ft>o < co[ < co' 2 < co' 3 < • • • < 00. 



(62) 



The transitions ejj D have the same sign for all even (odd) 
singular frequencies, which is opposite to that of odd (even) 
singular frequencies. 

The intersection of a singular line with kj = and fixed 

*° (0) = -ft)tan(ft)L+ ( Hft))) (63) 

discriminates between even and odd singular lines, too. 
For quasi-polynomials of retarded type (no infinity root 
boundaries), as ft) — >• 00, one group of the boundary straight 
lines diverges toward +°°, while the other towards — 00. Since 
for high frequencies, kp> and ep> share the same sign, the 
stable region must lie between even and odd boundary lines. 



Thus, starting from a sufficiently large singular frequency, 
the boundary lines become irrelevant. The same holds for 
the quasi-polynomials of neutral type. However, the infinity 
root boundaries ( |59| ) may impact the stable polygons: indeed 
these must lie between the two infinity root boundaries ( |6*0| . 
As a conclusion, stable inner polygons are determined by 
the low frequency boundaries and infinity root boundaries 
(if any). For instance, it was shown that for PID control 
of a first order proportional term G = K/(Ts+ \)e~ Ls , the 
relevant singular frequencies are the first two ones and the 
two infinity root boundaries. 

It is difficult to state a rule, which would precisely 
discriminate the relevant frequency range in the general case. 
However, some estimates are still thinkable. The relevant 
frequency range should comprise the singular frequencies, 
which correspond to the minima and maxima of the stable 
fcp-intervals. Another helpful rule of thumb is to discriminate 



the region, where kp(co), Fig. 1(a) oscillates with an almost 
fixed period as settled in ( |55) . 

C. Stable kp intervals 

In this section the solution of kp— problem is extended to 
time-delay systems. 

Theorem 5: Consider the quasi-polynomial ( |5Q| ). Assume 
that A (s) has J non-zero zeros on the imaginary axis, P right- 
hand-side zeros and a zero of order Jq at s = 0. // ( |50[ ) is 
Hurwitz-stable for a fixed kp = k' p , then a k € N exists, such 
that for I >k, I £ N, the number of singular frequencies Z in 
the interval < ft) < [2ln + 8)/L corresponding to kp = k' p 
is 

E(4l+N-M + 2P-J-l)+E(Jo) 



Z> 



(64) 



where 8 is chosen such that the principal term does not 
vanish at CO — {±2lK + 8)/L. 

For the proof of this theorem is refer to [Bajcinca (2004)]. 

Example 7. Consider PID control of the plant 

rM = -/-7.9 3 -2s +\ 005s 

U (s+l)(s + 2)( S + 3)(s + 4)(s 2 +s+l) ■ { ' 

Note that N = 7, M = 4 and P = 1 (since a right-hand-sided zero 
of A(s) exists at s = 0.3483). According to Theorem [5] we need 
to find a sufficiently large k, such that within any interval < 
co < 40/tt -f- 8, with / > k, at least £(4/ + 6)/2 = 2r + 3 singular 
frequencies are available. It can be easily checked that already for 
k = 1 and 8 = it the condition is fulfilled within — 24 < k P < 6.0693. 
Fig. [TO] shows the set of all PID stabilizers for the plant j65}. This 
example illustrates two interesting situations: first, for —3.7671 < 
kp < 4.6807 the stable region includes two separated polygons, and 
second, one of the polygons closes at a vertex. 

VIII. Conclusion 

Fast computational methods of the set of all PID con- 
trollers for linear continuous-time, discrete-time and time- 
delay systems are proposed in this article. The driving force 
of the theory is the fact that non-convex stability regions 
can be built up easily by convex polygonal slices. The high 
computational speed results due to inspection of conditions 
at a relatively low number of (singular) frequencies. Only 




Fig. 10. PID stability region, Example 7 
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the results of the control group at DLR are surveyed here. 
A software tool called ROBSIN originated on that basis. 

In author's opinion, the proposed design approach is espe- 
cially elegant in the discrete-time domain, and of particular 
interest for time-delay systems. A powerful feature is the fact 
that, in principle, in all cases all design lines apply also when 
simultaneous stabilization of a set of plants is considered. 
This paves the basis for robust design of PID controllers. 
Indeed, while in all derivations of the article just a single 
representant is assumed, the extensions to the situation with 
a finite number of representants is straightforward. When 
applying the £p-criterion, one would have to search for the 
intersection of stable &p— intervals of each representant. And, 
inner-polygon candidates should provide the simultaneous 
stability for all representants. Yet some important issues, 
particularly those involving time-delay systems, remain open. 
For instance, it is not clear how to discriminate the fre- 
quency range with relevant singular frequencies for quasi- 
polynomials. 
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